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Abstract PSR B0540—69 has braking index measurement in its persistent state: n = 

2.129 ± 0.012. Recently, it is reported to have spin-down state changes: a suddenly 36% 
increase in the spin-down rate. Combining the persistent state braking index measurement 
and different spin-down states, PSR B0540—69 is more powerful than intermittent pulsars 
in constraining pulsar spin-down models. The pulsar wind model is applied to explain the 
variable timing behavior of PSR B0540—69. The persistent state braking index of PSR 
B0540—69 is the combined effect of magnetic dipole radiation and particle wind. The 
particle density reflects the magnetospheric activity in real-time and may be responsible 
for the changing spin-down behavior. Corresponding to the 36% increase in the spin- 
down rate of PSR B0540—69, the relative increase in the particle density is 88% in the 
vacuum gap model. And the model calculated braking index in the new state is n = 1.79. 

Future braking index observation of PSR B0540—69 in the new spin-down state will be 
very powerful in distinguishing between different pulsar spin-down models and different 
particle acceleration models in the wind braking scenario. The variable timing behavior 
of PSR J1846—0258 is also understandable in the pulsar wind model. 
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1 INTRODUCTION 

PSR B0540—69, known as the “Crab Twin”, is a young radio pulsar with spin-down parameters 
v « 19.727 Hz, v ss —1.86 x 10 -10 Hzs -1 (Marshall et al. 2015) and braking index n = 2.129 ± 0.012 
(Ferdman et al. 2015). Its characteristic magnetic field is about 10 13 G at the magnetic polefj Only 
two glitches with relative small changes in spin-down parameters were reported (Zhang et al. 2001; 
Cusumano et al. 2003; Livingstone et al. 2005; Ferdman et al. 2015). Recently, a persistent and un¬ 
precedented increase in the spin-down rate of PSR B0540—69 was observed: the relative increase in the 
spin-down rate is 36% which is orders magnitude larger than the changes induced by glitches (Marshall 
et al. 2015). Another pulsar PSR J1846—0258 was also reported to have variable timing behaviors: a net 
decrease in the spin frequency (Au « —10 4 Hz) after the large glitch (Livingstone et al. 2010) and a 
lower braking indexn = 2.19 ± 0.03 (Livingstone et al. 2011; Archibald et al. 2015b)than its persistent 
state value n = 2.65 ± 0.01 (Livingstone et al. 2006). 

The spin-down behavior of pulsars can be described by the power law: 


v = -Cv n 


( 1 ) 
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Table 1 Comparison of spin-down parameters of PSR B1931 +24, PSR B0540—69, and PSR 
J1846 -0258. The intermittent pulsar PSR B1931+24 has different spin-down states without 
any braking index information at present. PSR B0540—69 has both different spin-down states 
and the persistent state braking index measurement. PSR J1846—0258 is reported to have a 
variation of braking index. 


Pulsar name 

iz(Hz) 

P(Hzs 

-1 

) 


braking index 

B1931+24 (off) a 

1.229 

-10.8 

X 

10" 

-15 

? 

B1931+24 (on) a 

1.229 

-16.3 

X 

10" 

-15 

? 

B0540-69 (low) b 

19.727 

-1.86 

X 

10" 

-10 

2.129 6 

B0540-69 (high) b 

19.701 

-2.53 

X 

10" 

-10 

? 

J1846—0258 (persistent state) d 

3.08 

-6.72 

X 

10" 

-11 

2.65 

J1846—0258 (after glitch) 6 

3.06 

-6.65 

X 

10' 

-11 

2.19 


(a) : From Kramer et al. (2006). The on state has larger spin-down rate than the off state. 

(b) : From Marshall et al. (2015). “Low” means the previous spin-down state and “high” means the new spin-down 
state with a higher spin-down rate. 

(c) Mean value of braking index (Feldman et al. 2015). 

(d) From Livingstone et al. (2006). 

(e) From Archibald et al. (2015b). 


where v and v are respectively the spin frequency and frequency derivative, C is usually taken as a 
constant and n is the braking index. The braking index is defined accordingly: 


where v is the second derivative of spin frequency. The braking index reflects the pulsar braking mecha¬ 
nism (Tong 2015). In the magneto-dipole braking model, a pulsar rotates uniformly in vacuum v oc 
The expected braking index is three which is not consistent with the observations (Lyne et al. 2015). Like 
the intermittent pulsar (Kramer et al. 2006), PSR B0540—69 also has two different spin-down states. 
For the intermittent pulsar PSR B1931+24, people tried to measure its braking index during the on and 
off state (Young et al. 2013). Now this aim has been partially fulfilled by PSR B0540—69 which has 
not only different spin down states but also braking index measurement for the persistent state (“low” 
spin-down rate state), see Table Q] Therefore, it can put more constraints on pulsar spin-down models. 
Any candidate model should explain both the braking index during the persistent state and the variable 
spin-down rate. 

Previously, the pulsar wind model (Xu & Qiao 2001) is employed to explain the spin-down behavior 
of intermittent pulsars (Li et al. 2014) and braking index of the Crab pulsar (Kou & Tong 2015). In the 
following, it is shown that both the persistent state braking index and varying spin-down rate of PSR 
B0540—69 are understandable in the wind braking model. The varying spin-down rate is due to a vari¬ 
able particle wind. And the varying braking index of PSR J1846—0258 is caused by a changing particle 
density. The pulsar wind model and calculations are listed in Section 2. Discussions and conclusions are 
presented in Section 3 and Section 4, respectively. 


2 VARIABLE TIMING BEHAVIOR OF PULSARS CAUSED BY A VARYING PARTICLE 
WIND 


2.1 Description of the pulsar wind model 


Pulsars are oblique rotators in general. The perpendicular and parallel magnetic dipole moment may 
respectively relate to the magnetic dipole radiation and particles acceleration (Xu & Qiao 2001; Kou & 
Tong 2015): 


E d 


2 fi 2 n 4 


sin 2 a. 


3c 3 


(3) 
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2 ^ 2 Q 4 
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o 2 

ok-— cos a, 
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(4) 


where p = 1/2 BR 3 is the magnetic dipole moment (B is the polar magnetic field and R is the neutron 
star radius), c is the speed of light, and a is the angle between the rotational axis and the magnetic 
axis (i. e., inclination angle), fi = 2-kv is the angular velocity of the pulsar, r p = R(RQ/ c ) 1 / 2 is the 
polar cap radius, p e = Kpc.j is the primary particle density where poj = HB / (2-c) is the Goldreich- 
Julian charge density (Goldreich & Julian 1969) and k is the dimensionless particle density, A</ is the 
corresponding acceleration potential of the acceleration region, and A<f> = //0 2 /c 2 is the maximum 
acceleration potential for a rotating dipole (Ruderman & Sutherland 1975). The pulsar rotational energy 
is consumed by the combined effect of magnetic dipole radiation and particle acceleration (Xu & Qiao 
2001 ) 

2 u 2 n 4 

- mn = (5) 

where I = 10 45 g cm 2 is the moment of inertia, and 


77 = sin 2 a + 3nA(f>/ A<1> cos 2 a. 


( 6 ) 


The spin-down behavior can be expressed as: 


fi = - 


2/i 2 fl 3 


(7) 


According the equation ©, the braking index in the pulsar wind model can be written as (Xu & 
Qiao 2001): 


ft dr] 

n — 3 4 ~77v 

77 ail 


( 8 ) 


The exact expression of 77 (equation ©I depends on the particle acceleration potential. The vacuum gap 
model (Ruderman & Sutherland 1975) is taken as an example to show the calculation process and 


rj = sin 2 a + 4.96 x 10 2 «.B 12 8 ' / 7 f2 15 ^ 7 cos 2 a, (9) 

where B 12 is the magnetic field in units of 10 12 G (Kou & Tong 2015). For other acceleration models, 
the corresponding expressions of 77 are listed in Table 2 in Kou & Tong (2015). 


2.2 On the variable timing behavior of PSR B0540—69 

A generic picture for the variable timing behavior of PSR B0540—69 and PSR J1846—0258 is: a glitch 
may occurred during the observations, as that in PSR J1846—0258 (Livingstone et al. 2010). This small 
glitch may have been missed in the case of PSR B0540—69. This glitch may induce some magneto- 
spheric activities, e.g., outburst (Gavriil et al. 2008). The particle outflow will be stronger during this 
process. This will cause the pulsar to have a larger spin-down rate (Marshall et al. 2015). After some 
time, a larger spin-down state will result in a net spin-down of the pulsar compared with previous timing 
solutions (Livingstone et al. 2010). The braking index will be smaller since the particle wind is stronger 
(Wang et al. 2012a). When the pulsar magnetosphere relax to its persistent state, if the particle density 
is still varying with time k = n (i), the braking index will be different with the persistent state while 
the change in spin-down rate may be negelected (Livingstone et al. 2011; Archibald et al. 2015b; Kou 
& Tong 2015). From previous observations of PSR J1846—0258, its pulse profile has no significant 
variations before, during and after the outburst (Livingstone et al. 2010, 2011; Archibald et al. 2015b). 
The observations of magnetar IE 1048.1 —5937 showed that the pulsed flux may not be a good indicator 
of the magnetospheric activities (Archibald et al. 2015a). The variation of total X-ray flux is needed. 
Therefore, the enhanced spin-down rate in PSR B0540—69 without changes in pulse profile and pulsed 
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Table 2 Braking indices of PSR B0540—69 in the new state in all the acceleration models. 


Acceleration models 

VG(CR) 

VG(ICS) 

SCLF(II,CR) 

SCLF(I) 

OG 

CAP 

NTVG(CR) 

NTVG(ICS) 

Braking index 

1.79 

1.87 

1.90 

1.79 

1.38 

1.83 

1.90 

1.86 


Notes: See Table 2 of Kou & Tong 2015 for the meanings of the acceleration models abbreviations. The minimum 
braking index of the SCLF (II ICS) is 2.4 (Li et al. 2014) which is larger than the persistent braking index of PSR 
B0540—69: n = 2.129. It means that the SCLF (II, ICS) model can be ruled out or can not exist alone to 
accelerate particle in the magnetosphere of PSR B0540—69 (Wu et al. 2003; Li et al. 2014). 


flux is not unusual (Marshall et al. 2015). The reason may be that the geometry of the pulsar is un¬ 
changed during the magnetospheric activities. This may result in a constant pulse profile and pulsed 
flux. 

For PSR B0540—69, giving the persistent state spin-down parameters v = 19.727 Hz, v = —1.86 x 
10 -10 Hz s -1 (Marshall et al. 2015), and inclination angle a = 50° (the best fitted value given by Zhang 
& Cheng 2000). Parameters of magnetic field B = 10 13 G and re = 834 can be calculated (by solving 
equations (0 and ©) corresponding to the observed braking index n = 2.129 ± 0.012 (Ferdman et al. 
2015). The calculated re = 834 means that the particle density is 834 times the Goldreich-Julian charge 
density which is consistent with previous conclusions (Kou & Tong 2015 and references therein). 

The spin-down rate of PSR B0540—69 has increased by 36% in the new spin-down state (Marshall 
et al. 2015). In the pulsar wind model, the variation of the spin-down rate is caused by a different particle 
density 


H ?y(re ) 

Cl r]{n) ’ 


GO) 


where Cl and r](n ) correspond to the new spin-down state. A larger particle density will result in 
a higher spin-down rate (equation® and ®). Figure Q] and figure [ 2 ] shows respectively the normal¬ 
ized spin-down rate v /re and braking index as function of normalized particle density re /re for PSR 
B0540—69 in the vacuum gap model. As shown in figure® the spin-down rate increases as the particle 
density increases. An increase in the particle density of 88 % will result in the 36% increase in the spin- 
down rate. As particle density increase, the braking index will decrease because the effect of particle 
wind component is increasing (figure®. When the particle density increases to 1.88 times the previ¬ 
ous value, the braking index decreases to 1.79, the relative change is 15.7%. And the corresponding 
frequency second derivative will be re = 5.83 x 10 -21 Hz s -2 (equation ®). 

Calculations in all the acceleration models are also made. The same conclusion is obtained 
from these models: an increasing particle density results in the increase in spin-down rate. For PSR 
B0540—69, corresponding to the observational 36% relative increase in spin-down rate, the relative 
increase in the particle density in all these models ranges from 72% to 154%. The second frequency 
derivative ranges from 4.5 x 10 " 21 Hz s -2 to 6.15 x 10 ~ 21 Hz s~ 2 . Braking indices in the new state in 
all these acceleration models are list in Table ® If the conversion efficiency of particle energy to X-ray 
luminosity is unchanged (Becker 2009), the total X-ray luminosity may also have increased by the same 
factor. 


2.3 On the variable timing behaviors of PSR J1846 0258 

Spin-down parameters and persistent state braking index of PSR J1846—0258 are respectively: re = 
3.08 Hz and re = —6.72 x 10~ 11 Hzs _1 and n = 2.65 ± 0.01 (Tabl® (Livingstone et al. 2006). 
In the pulsar wind model, corresponding to the observational braking index, the magnetic field B = 
1.25 x 10 14 G and particle density re = 28 are calculated in the vacuum gap model with an inclination 
angle 45° (a inclination angle of 45° is chosen in the following calculation^). Such a magnetic field 
is comparable with the characteristic magnetic field 9.7 x 10 13 G at the poles and much larger than 


2 


There is no observational or best fitted inclination angle given. 
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Fig. 1 The normalized spin-down rate as function of the normalized particle density for PSR 
B0540—69 in the vacuum gap model. The dashed line is the spin-down rate in the persistent 
state. The dotted line is the new spin-down rate which is 1.36 times the persistent state spin- 
down rate (Marshall et al. 2015). 



normalized particle density 


Fig. 2 Braking index as function of the normalized particle density for PSR B0540—69 in the 
vacuum gap model. The dashed line is the persistent state braking index 2.13 (Ferdman et al. 
2015). The dotted line is the braking index 1.79 which is predicted by the increased particle 
density. 


magnetic fields of normal pulsars. Then it is not surprising that magnetar activities can be observed in 
this source (Gavriil et al. 2008). 

Variable timing behavior of a net decrease in the spin-down frequency (A v « -10' 1 Hz) was 

detected for PSR J1846—0258 after a larger glitch (Livingstone et al. 2010). The correspondingly rel¬ 
ative increase in the spin-down rate is about 7% (Az> = —4.82 x 10 -12 Hzs -1 during an epoch 240 
days when phase coherency is lost). Such an increase in spin down rate may also be caused by a larger 
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particle density. Just like the calculation for PSR B0540—69, in the vacuum gap model of the pulsar 
wind model, a 44% increase in the particle density results in the 7% increase in the spin-down rate. 
The braking index will also be smaller during this enhanced spin-down epoch. However, braking index 
measurement is only available long after the glitch when the timing noise is greatly reduced. 



Fig. 3 The braking index of the PSR J 1846—0258 as function of r K in the vacuum gap model. 
The dashed line is the persistent state braking index 2.65 (Livingstone et al. 2006). The dotted 
line is the smaller braking index 2.19 measured after the glitch (Archibald et al. 2015b). 


A lower braking index 2.19 ± 0.03 is detected after the glitch (Livingstone et al. 2011; Archibald 
et al. 2015b) which is significantly smaller than its persistent state value n = 2.65 ± 0.01 (Livingstone 
et al. 2006). In the pulsar wind model, it can be understood by a time varying particle density k = n{t) 
(similar to the Crab pulsar, Kou & Tong 2015): 


n dr] 

n = 3 + “377 
r] ail 


k dr/ t c 
r] dur ^ 1 


( 11 ) 


where r c = --=4 is the characteristic age, r K = is the typical variation timescale of the particle 

density. An increasing particle density (r K > 0 or k > 0) lead to a smaller braking index. From figure 
[3 the braking index is insensitive to t k when it is much larger that 1000 yr, but decrease sharply when 
r K is comparable with the characteristic age (about 700 yr). The changing rate of the particle density 
k = 1.68 x 10 -9 s -1 will result in a braking index of 2.19. During the epoch from MJD 55369 to MJD 
56651, the particle density has increased 0.66%. The relative increase in spin-down rate is 0.1% which 
is very small. Therefore, the changing particle density will mainly result in a different braking index 
while not affecting the spin-down rate. This is the difference between the variable timing behavior of 
PSR 0540-69 and PSR J1846-0258. 


3 DISCUSSIONS 

Observations of intermittent pulsars (Kramer et al. 2006) and measurement of braking indices (Lyne 
et al. 2015) help to distinguish between different pulsar spin-down mechanisms. The variable spin- 
down rate of PSR B0540—69 combined with its persistent state braking index measurement is more 
powerful than intermittent pulsars in constraining different models. In the magneto-dipole radiation 
model v oc /u 2 sin 2 a/1 x i/’\ in order to explain the braking indices (n < 3) of eight young pulsars, 
an increasing inclination angle (Lyne et al. 2013), increasing magnetic field (Espinoza et al. 2011) or 
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decreasing moment of inertia (Yue et al. 2007) is expected. The variable spin-down rate may be induced 
by the change of inclination angle, magnetic field or moment of inertia. Corresponding to the 36% 
relative increase in the spin-down rate, the relative change should be: 26% increase in inclination angle, 
or 17% increase in magnetic field, or 26% decrease in moment of inertia. It seems impossible to achieve 
such huge changes during a short timescale (about 14 days, Marshall et al. 2015). A change in inclination 
angle is unlikely since the pulse profile did not change significantly (Marshall et al. 2015). An increase 
in magnetic field will require an increase of magnetic energy by 36%, about 10 42 erg. It is unlikely there 
is such a huge amount of energy injection. A decrease of moment of inertia will require a decrease of 
neutron star radius. During this process, a huge amount of gravitational energy will be released (Zhou 
et al. 2014), about 10 52 erg, which is again unlikely. 

Previous models for the spin-down behavior of intermittent pulsars (Beskin & Nokhrina 2007; Li et 
al. 2012) may also be applied to the variable spin-down rate of PSR B0540—69. However, the expected 
braking index is three in Beskin & Nokhrina (2007) and the magnetohydrodynamical simulations (Li 
et al. 2012). Considering the effect of pulsar death or evolution of inclination angle, the braking index 
will be larger than three (Contopoulos & Spitkovsky 2006; Philippov et al. 2014). Therefore, these 
models should be modified before they can explain both the persistent state braking index and variable 
spin-down rate of PSR B0540—69. 

There are several models designed for magnetar spin-down which may also be employed to the 
case of PSR B0540—69. The magnetar spin-down may be dominated by a particle wind (Harding et 
al. 1999). The calculations in Harding et al. (1999) is equivalent to assuming each outflow particle can 
attain the maximum acceleration potential of a rotating dipole (Tong et al. 2013). This wind braking 
model of magnetars was employed by Kramer et al. (2006) to explain the spin-down behavior of the 
first intermittent pulsar PSR B1931+24. An additional particle outflow in the on state will result in a 
larger spin-down rate. The rotational energy loss is related to the particle wind luminosity L p as oc s/L^, 
(Harding et al. 1999). A particle wind luminosity 85% larger will result in a spin-down rate 36% larger. 
The particle wind luminosity is related to the polar cap radius R pc and magnetospheric opening radius 
r op en as L p oc R pc oc r~ pen (Harding et al. 1999). Therefore, the magnetospheric opening radius will be 
26% smaller. However, there are several problems when applying the wind braking model of magnetars 
to the case of normal pulsars: (1) In the wind braking model of magnetars, a strong particle wind is 
assumed. The effect of magnetic dipole radiation is neglected. This may be appliable to the case of 
magnetars whose emissions are dominated by magnetic energy output (Tong et al. 2013). However, in 
the case of normal pulsars (including intermittent pulsars) the effect of dipole radiation may not be 
neglected. (2) In the case of strong particle wind, the braking index is n = 1 (Tong et al. 2013). This 
is not consistent with the braking index of pulsars (Lyne et al. 2015). (3) When applying to the case of 
intermittent pulsars (Kramer et al. 2006; Young et al. 2013), pure magnetic dipole braking is assumed 
for the off state. This may be valid for the case of intermittent pulsars whose radio emissions are stopped 
in the off state (Li et al. 2014). However, this assumption can not be applied to the persistent spin-down 
state of PSR B0540—69 which still has multiwave emissions. 

The twisted magnetosphere model of magnetars (Thompson et al. 2002) showed that the effective 
magnetic field will be larger for a larger twist. If the magnetosphere of PSR B0540—69 is twisted by 
a glitch, then it will also result in a larger spin-down rate. However, the twisted magnetosphere will 
relax back to the pure magnetic dipole case in several years (Beloborodov 2009). During this process, 
the neutron star X-ray luminosity, spin-down rate will both decrease with time. For PSR B0540—69, its 
high spin-down state has lasted more than 3 years (Marshall et al. 2015). This is inconsistent with the 
expectation of the twisted magnetosphere model. 

There are also external models for the braking index or intermittent pulsar spin-down behavior, e.g., 
the fallback disk model (Liu et al. 2014 and references therein; Li et al. 2006). However, these external 
models are hard to verify or falsify. Furthermore, accretion will halt the magnetospheric activities. In 
the presence of accretion, it may be difficult to reconcile with the radio emissions in PSR B0540—69 
and in other pulsars with braking index measured. 

Observations of the timing behavior and pulse profile of some pulsars indicate that the v modula¬ 
tion and the pulse-shape variation are correlated, e.g., PSR B0910+16 (Perera et al. 2015a) and PSR 
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B1859+07 (Perera et al. 2015b). The connection indicates that both these phenomena are magneto- 
sphereic origin (Lyne et al. 2010). Theory of a variable particle density in the magnetosphere is success¬ 
fully applied to explain the spin-down behavior and emission property of the intermittent pulsar (Kramer 
et al. 2006; Li et al. 2014). The mode changing and nulling pulsar may be understood similarity because 
of the detection of variation in spin-down rate of PSR .11717 —4054 (Young et al. 2015) and a weak 
emission state in addition to its bright and nulling states of PSRs J1853+0505 and J1107-5907 (Young 
et al. 2014). For PSR B0540—69, the increase of particle density in the magnetosphere will change the 
spin-down rate, as well as the pulse profile. The radio giant pulse of PSR B 0540—69 (Johnston et al. 
2004) may be caused by a larger out-flowing particle density. Besides, different emission models (core, 
cone and patch) are applied to explain the variable mean pulse profiles (Lyne & Manchester 1988). And, 
it is predicted that the corotation of magnetosphere with pulsar may also affects the emission properties 
(Wang et al. 2012b). Hence, the nonuniform distribution of particles between these core and conal com- 
ponnets will change the pulse shape also. For PSR B0540—69, the pulse profile has a broad double peak 
which can be described with two Gaussians with a phase separation of 20% (de Plaa et al. 1993).We 
could emphasize that: (i) if the particles distribute uniform, the increase in out-flowing particle density 
may result in the increase of pulse intensity, the ratio of these two component keeps constant; (ii) if 
the particles distribute nonuniform, both the pulse intensity and the ratio will change; (iii) the coherent 
manner of radiated particles may affect the pulse shape as well. The pulse profile and total flux in the 
high spin-down state are needed to compare with them in the low spin-down state. 

4 CONCLUSIONS 

The pulsar wind model is applied to explain the variable timing behavior of PSR B0540—69 and PSR 
J1846—0258. Both the persistent state braking index and the variable spin-down rate of PSR B0540—69 
are understandable. A larger particle density will result in an increase in the spin down rate and predicts 
a smaller braking index. And an increasing particle density will lead to a lower braking index. For 
PSR B0540—69, in the vacuum gap model, corresponding to the 36% increase in the spin-down rate, 
the relative increase in particle density is 88%. And the braking index decreases to 1.79. The same 
conclusion is obtained for the different acceleration models. Since it has both a variable spin-down rate 
and persistent state braking index measured, PSR B0540—69 is very powerful in constraining different 
pulsar spin-down mechanisms. Future observation of braking index in the new spin-down state will 
provide further test on different spin-down models and different particle acceleration models in the 
wind braking scenario. For PSR J1846—0258, the variable timing behavior of a net decreasing in spin 
down frequency (Az/ ss —10 -4 Hz) can be understood similarly. And a changing rate of particle density 
k = 1.68 x 10 -9 s -1 will result in the lower barking index 2.19. 
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